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ABSTRACT 

Simulations of the formation of late-type spiral galaxies in a cold dark matter (ACDM) universe have 
traditionally failed to yield realistic candidates. Here we report a new cosmological A^-body/smooth 
particle hydrodynamic (SPH) simulation of extreme dynamic range in which a close analog of a Milky 
Way disk galaxy arises naturally. Termed "Eris" , the simulation follows the assembly of a galaxy halo 
of mass Mvir = 7.9 x 10^^ M0 with a total of = 18.6 million particles (gas -I- dark matter -I- stars) 
within the final virial radius, and a force resolution of 120 pc. It includes radiative cooling, heating 
from a cosmic UV field and supernova explosions (blastwave feedback), a star formation recipe based 
on a high gas density threshold (nsF = 5 atoms cm~'^ rather than the canonical nsF =0.1 atoms 
cm~^), and neglects any feedback from an active galactic nucleus. Artificial images are generated to 
correctly compare simulations with observations. At the present epoch, the simulated galaxy has an 
extended rotationally-supported disk with a radial scale length Rd = 2.5 kpc, a gently falling rotation 
curve with circular velocity at 2.2 disk scale lengths of V2.2 = 214 kms~^, an i-band bulge-to-disk 
ratio B/D = 0.35, and a baryonic mass fraction within the virial radius that is 30% below the cosmic 
value. The disk is thin, has a typical H i-to-stellar mass ratio, is forming stars in the region of the 
EgFR-Sui plane occupied by spiral galaxies, and falls on the photometric TuUy-Fisher and the stellar 
mass-halo virial mass relations. Hot (T > 3 x 10^ K), X-ray luminous halo gas makes only 26% of 
the universal baryon fraction and follows a "fiattened" density profile oc r~^-^^ out to r = 100 kpc. 
Eris appears then to be the first cosmological hydrodynamic simulation in which the galaxy structural 
properties, the mass budget in the various components, and the scaling relations between mass and 
luminosity are all consistent with a host of observational constraints. A twin simulation with a low 
star formation density threshold results in a galaxy with a more massive bulge and a much steeper 
rotation curve, as in previously published work. A high star formation threshold appears therefore key 
in obtaining realistic late-type galaxies, as it enables the development of an inhomogeneous interstellar 
medium where star formation and heating by supernovae occur in a clustered fashion. The resulting 
outflows at high redshifts reduce the baryonic content of galaxies and preferentially remove low angular 
momentum gas, decreasing the mass of the bulge component. Simulations of even higher resolution 
that follow the assembly of galaxies with different merger histories shall be used to verify our results. 
Subject headings: galaxies: evolution - halos - kinematics and dynamics - method: numerical 



1. INTRODUCTION 

The formation of realistic late-type spirals has been a 
long standing problem of galaxy formation in a ACDM 
universe. Within this framework, baryons condense at 
the center of dark matter halos and acquire angular mo- 
mentum through tidal to rques from nearby structures 
(Fall & Efstathiou 1980). A centrifugally-supported 
baryonic disk forms, with a size that depends on the 
fraction of the original angular momentum that is re- 
tained during the contraction. In numerical simulations 
of this process, however, a fundamental "angular mo- 
mentum problem" arises, as galaxies are produced with 
a baryonic component that is quite de ficient in angular 
mo mentu m compared to real spi rals (Navarro & Benz 
p9l| |N avarro & Steinmetz||2000|). Aside from artificial 
losses of angular momentum cause d by insufficient reso- 
lution and other numerical effects ( Okamoto et aH|2003 



Governato et al.|2004 Kaufmann et al.|2007[ ), this failure 
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has traditionally been traced back to the very nature of 
the hierarchical buildup of structures: dynamical friction 
transfers the orbital angular momentum of merging sub- 
structures to the outer halo, and causes the associated 
cold baryons to sink to the center of the proto-gala xy and 
form a spheroid rather than a disk (e.g. jMaller fc Dekel 



20021 



S'^popular solution to the angular momentum prob- 
lem envisions energy injection from supernovae (SNe) 
and evolving stars as a mechanism to prevent efficient 
gas cooling and condensation and to remove low angu- 
lar momentum material from the central part of galax- 
ies. Modern simulations with improve d resolution and 
more effective recipes for SN feedback (Robertson et al. 



2004^ 'Governato e t al.| [20071 JSc annap eco ai I [i^OIT^ 
Stins on et a l. 2010; |Piontek &: St cinmctz 2011| |Brooks 



et al.| |2011| niave yielded rotationally-supported disks 
with realistic exponential scale lengths, not only in galax- 
ies formed in relative isolation but also in those that 
are accreted by massive gro ups with a dominant cen- 
tral elliptical (Feldmann et al. 2010a b). They have also 
modified the standard picture of gas accretion and cool- 
ing onto galaxy disks: for galaxies up to Milky Way 
masses, gas acquired through filamentary "cold fiows" 
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TABLE 1 

Properties of the simulated galaxy 



Galaxy 


Mvir 




M, 


fb 


/"cold 


muM 


mspH 












B/D 


Rd 


Eris (z = 0) 


7.9 


238 


3.9 


0.121 


0.12 


9.8 


2 


120 


0.1 


7.0 


3.0 


8.6 


0.35 


2.5 


Eris (z = 1) 


5.4 


237 


2.9 


0.126 


0.40 


9.8 


2 


120 


0.1 


4.8 


2.0 


6.2 


0.30 


1.8 


ErisLT (z = 1) 


5.5 


308 


3.4 


0.158 


0.18 


9.8 


2 


120 


0.05 


4.9 


2.9 


8.3 


0.42 


1.4 



Note. — Columns 2, 3, 4, 5, and 6 list the virial mass (in units of 10^^ Mq), peak circular velocity (in kms~^), 
total stellar mass of the halo (in units of 10^" ^'^0)1 baryonic mass fraction, and cold (T < 3 X 10** K) gas fraction. 
Columns 7 and 8 list the mass resolution of individual dark matter and SPH particles (in units of 10* Mq), and 
columns 9 and 10 the spline gravitational force softening (in pc) and the star formation efficiency. Columns 11, 12, 
and 13 list the total number (in units of 10*') of dark matter, gas, and star particles within the virial radius of the halo. 
Columns 14 and 15 list the bulge-to-disk ratio and disk scale length (in kpc) estimated from the i-band photometric 
decomposition. 



that was never shock-heated to the halo virial tempera- 



all times (Brooks et al. 12009 


Keres et al.||2009 


Ceverino 


et al. 


201Up. Yet, these simulations typically continue 



curves that rise steeply towards the center: simulated 
disk galaxies fall exclusively in the SO or Sa category, 
leaving late-type spirals with negligible bulges and large 
disks with flat rotation curves - such as our own Milky 
Way - as an unsolve d puzzle. Two recen t exc eptions 
are the simulations of Agertz et al. (2011) and 'Gover-J 
nato et al. (2010). In the hrst, replicas of Sb/Sc galax- 
les with moderate bulges were obtained with a low effi- 
ciency of star formation that may implicitly mimic the 
bottleneck of the conversion of atomic gas into molec- 
ular, at the expense of producing stellar disks that are 
muc h more massive than expected at a given halo mass 
(e.g. Guo et al.|2010 |. In the second, a reahstic bulgeless 
dwarf galaxy with a shallow central dark-matter profile 
was generated by resolving the inhomogenous interstellar 
medium (ISM) and the process of energy injection from 
multiple SNe in clustered star forming regions. In this 
paper we extend the latter approach to massive galaxy 
scales, and present initial results from a new SPH cosmo- 
logical simulation of high dynamic range that includes 
radiative cooling, heating from a cosmic UV field, SN 
feedback, and a star form ation recipe based on a high 



gas density threshold as in Governato et al 
this last feature, we argue, that is key to t 



( |2010| ). It is 

le formation 



of a more realistic massive late-type spiral in ACDM. 

2. SIMULATION SETUP 

Dubbed "Eris" , the simulation described in this paper 
is part of a campaign of extreme resolution si mulations 
of the formation of Milky Way-sized galaxies (Diemand 
et al.|[200 7, 2008). It was performed m a Wilkinson Mi- 
crowave Anisotropy Probe 3-year cosmology, ^Im = 0.24, 
flA = 1- ^M, = 0.042, Ho = 73kms-iMpc-\ 
n — 0.96, as — 0.76, running the parallel, spatia lly and 
temporally ad aptive, treeSPH-code GASOLINE ( [Wads- 
ley et al.|2004 ) for 1.5 milion cpu hours. The target halo 
was identified at z = in a low-resolution, dark matter- 
only, periodic box of 90 Mpc on a side. It was choosen to 
have a similar mass as the Milky Way and a rather quiet 
late merging history, i.e. to have had no major mergers 
(defined as mass ratio > 1/10) after z = 3. New ini- 
tial conditions were then generated with improved mass 
resolution, centered around a Lagrangian sub-region of 



1 Mpc on a side, using the standard "zoom-in" tech- 
nique to add small-scale perturbations. High-resolution 
particles were further split into 13 million dark matter 
particles and an equal number of gas particles, for a fi- 
nal dark and gas particle mass of muM = 9.8 x 10'' M0 
and mspH — 2x 10^ M©, respectively. The gravitational 
softening length, eg, was fixed to 120 physical pc for all 
particle species from z = 9 to the present, and evolved 
as 1/(1 -I- z) from z = 9 to the starting redshift of z = 90. 

The version of the code used in this study in- 
cludes Compton cooling, atomic cooling, and metallicity- 
dependent radiative co oling at low temperatures 



(Mashchenko et al. 2006). A uniform UV background 
modihes the ionization and excitation state of t he gas and 
is implemente d using a modified version of the |Haardt &[ 
Madau ( 1996|) spectrum. Three parameters characterize 
the star formation and feedback recipes: (a) the star for- 
mation threshold tt-sf, (b) the star formation efficiency 
esF, and (c) the fraction of SN energy that couples to the 
ISM esN- Star formation occurs when cold (T < 3 x lO"' 
K), virialized gas reaches a threshold density hsf = 5 
atoms cm~^ and is part of a converging flow. It proceeds 
at a rate 



1.5 

gas 



(1) 



(i.e. locally enforcing a Schmidt law), where p* and pgas 
are the stellar and gas densities, and tdyn is the local 
dynamical time. We choose cqf = 0.1. An addi tional 
identical run with e sF = 0.05, the value adop ted in |Gov-| 
ernato et aT] ( |2007[ ) and [Brook et al] ( |2010[ ), yielded a 
galaxy with nearly identical structural properties, and 
will be discussed in a forthcoming paper. Each star 
particle is created stochastically with an initial mass 
= 6 X 10'^ Mq, and the gas particle that spawns the 
new star has its own mass reduced accordingly. A star 
particle represents a simpl e stellar population with its 
own age, metallicity, and a Kroupa 32001 ) initial stellar 
mass function (IMF) . Each SN deposits metals and a net 
energy of esN x 10^^ ergs into the nearest neighbor gas 
particles, with esN = 0.8 (the same value adopted in pre- 
vious simulations). The heated gas has its cooling shut 
off until the end of the snowplow phase of the SN blast- 
wave, which is set by the local gas density and tem pera- 



ture and by th e total amount of energy injected E (Stin- 



old touu' 
Re 



son et al^2006[). For the typical ISM conditions at thresh 



Tnthis study, this translates into regions of size 
pc heated by individual SNe and having 
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their cooling shut off for a timescale is ^ 5 x 10^ E^^^ yr, 
where E^i = ergs. The energy injected by many 

SNe adds up to create larger hot bubbles and longer shut- 
off times. No feedback from an active galactic nucleus 
was included. 

The adoption of a density threshold for star formation 
that is 50 times higher than in many previous lower- 
resolution studies is possible owing to the high mass and 
spatial resolution of this run, which resolves the giant 
cloud complexes where star formation actually occurs in 
the ISM and the true scale height of the neutral atomic 
ISM. In particular, the local Jeans length corresponding 
to our density threshold (for T = 10^ K, a lower bound 
on the typical temperature of the cold gas in the sim- 
ulations) is resolved with more than 5 SPH smoothing 
lengths, thus pre venting artificial fragmentation ( .Bate| 
& Burkert 19971. While not as high as the value of 
nsF = 100 atoms cm~'^ used in Governato et al. (20101 
dwarf galaxy simulation, whose particle mass was sub- 
stantially lower and allowed to resolve the Jeans length 
of star forming gas at much higher densities, our star 
formation threshold is still large enough to allow the de- 
velopment of a clumpy, inhomogeneous ISM with more 
localized energy injection by multiple overlapping SN ex- 
plosions. This allows galactic pressure-driven outflows 
to develop and remove low-angular momentum material. 
To demonstrate the important role of the star formation 
threshold on the structural properties of massive galax- 
ies, we have run a low-threshold twin simulation (termed 
"ErisLT") with ngp = 0.1 atoms cm~^. We have kept all 
the other simulation parameters fixed (same mass and 
spatial resolution and identical feedback scheme) except 
for the star formation efficiency parameter, esF, which 
was lowered from 0.1 (Eris) to 0.05 (ErisLT) to match 
the observed normal ization of the star form ation density 
in local galaxies (see Governato et al.l2010[ ). ErisLT was 
stopped at redshift 0.7 in order to limit the computa- 
tional burden. 

3. RESULTS 

In this first paper we focus on the final z = state 
and properties of the galaxy and on the comparison with 
observational constraints. The main characteristics of 
the simulated galaxy are listed in Table 1. 

3.1. Structural Parameters 

At the present epoch, Eris is a massive, barred, late- 
type spiral with structural properties consistent with 
those of Sb/Sbc galaxies, of which our own Milky Way 
is an example. It has a virial radius of i?vir = 239 
kpc (defin ed as the radius enclos ing a mean density 



of 93/9crit, Bryan & Norman 1998), total mass M, 



7.9 X 10" M (7), spin parameter A = 0.019 (a la [Bullock 
et al. 2001), and 7.0, 3.0, and 8.6 million dark matter, 
gas, and star particles within i?vir, respectively. The 
minimum smoothing length for gas particles is 5 times 
smaller than the force softening. The total mass en- 
closed within 60 kpc is Af<6o = 3.3 x 10" Mq. The 
rotation curve, shown in the left panel of Figure [ij has 
a peak circular velocity of Vpeak = 238 kms~^ (reached 
at 1.34 kpc) and a value at 8 kpc (the solar circle) of 
Vc^Q = 206 kms~^. Its overall shape out to 20 kpc, 
including the peak in the central bulge-dominated kpc. 



is reminiscent of the r ecent reconst ruction of the Milky 
Way rotation curve by jSofue et al.| p009). The circular 
velocity decreases gently to distances of 60 kpc from its 
value at the solar radius, in agreement with observations 
of blue hori zontal-branch halo stars in the Sloan Digital 
Sky Survey ( |Xue et al.|2008) . The measured K,©, A/<go, 
and Mvir agree wittiin the errors with the values of 
Vc,0 = 221 ± 18 kms-i, M<6o = 4.0 ± 0.7 x 10" M©, 
and Mvir = ^-Q^q.I x 10^^ y^^ derived recently for the 
Milk y Way using the narr ow GD-1 stream of stars (for 



Vr 



c,0> 



tracers (tor M^qq and M 



Koposov et al. 2010|) and halo stars a s kinematic 
]Xue et al.|[2008 ). 



3.2. Brightness Profile 

To correctly compare simulations with observations 
we created artificial images of our simulations and from 
them measured photometric bulge-to-disk ratios and disk 
scale lengths. The mock images were created u sing the 
radiation transfer code Sunrise (Jonsson 2006), which 
produces spectral energy distributions using the age and 
metallicities of each simulated star particle, and takes 
into account the three-dimensional effect of dust repro- 
cessing. The results for a Kroupa IMF are shown in Fig- 
ure [2j A 2D photometric decomposition was performed 
on the dust-redden ed i-band light dis tribution with the 
Galfit program (Peng et al. 2002). At the present 
epoch, the total z-band magnitude is Mi = —21.7, and a 
stellar disk with a scale length Rd = 2.5 kpc dominates 
the light distribution (Fig. |3|. The disk scale length is 
comparable to the value Rd — 2.3 ± 0. 6 kpc, adopted 
for the Milky Way in the compilation by [Hammer et ah] 
( [2007[) , and with the scale length of the Milky Way thin 
disk, 2.6 kpc, as traced by M dwarfs in the solar neigh- 
borhood (Juric et al. 2008). Its value also agrees with 



the s caling relations of spiral galaxies ( Courteau et al 
2007|). The SDSS u - .g = 1.03 mag and g - r = 0.49' 



mag integrated colors, obtained directly from the Sun- 
rise images, fall within la of the mean optical colors 



of lat e-type galaxies as luminous as Eris (Blanton et al. 
2003|. The "downbending" observed in Eris' brightness 



exponential profile at about 4 disk sca le lengths appears 



to be characteristic of late-type spirals (Pohlen & Trujillo 
20061). As in the sample of trun cated late-type spirals of 



akos, Tr ujillo, k Pohlen[ ([20081 , there is no break in the 
stellar surtace mass density profile of Eris: rather, Eris' 
stellar age profile shows a "U shape" with a minimum 
of 6 Gyr at the break radius, explaining the origin of 
the break as a radial change in stellar population likely 
caused by the stochastic radial migration of young stars 
from the inner parts of the disk to the outskirts (Roskar] 



et al. 



2008) 



Eris bulge-to-disk ratio (as determined by a two- 
component fit to the i-band surface brightness pro- 
file), B/D = 0.35, is also typical of Sb spirals, which 
are characterized by a median (±68/2 per cent) value 
\ogB/D = — 0.53to'3o: and of many Sbc galaxies, which 



have log B/D = -OMtllt ([Graham & Worley 2008[). A 



three-component decomposition (disk-|-bar-|-bulge) will 
lower the B/D ratio further. The bulge Sersic index. 
Us = 1.4, is indicative of a "pseudobulge" rather than a 
classical one: according to Weinzirl et al. ( |2009 ), ~ 3/4 
of all bright spirals have low < 2 bulges. Eris large fi- 
nal disk (disk-to-total ratio D/T = 0.74) is not typically 
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Fig. 1. — Left panel: The rotation curve of the simulated Milky Way-sized galaxy ("Eris") at 2 = 0. The figure shows the contributions 
to the circular velocity Vc = GM{< r)/r of the various mass components: dark matter {long-dashed curve), stars (short-dashed curve), 
gas (dot-short dashed curve), and total (solid curve). The data points show two realizations of the rotation curve of the Milky Way from 
observations of blue horizontal-branch halo stars in the Sloan Digital Sky Survey (Xuc et al.|2008| ), and have been offset slightly from each 
other in radius for clarity. Right panel: The total inner rotation curve at z = 1 for the fiducial mgh-threshold simulation (Eris, solid line) 
and for the low-threshold (ErisLT, dot-short dashed line) twin run. The short-dashed line shows Eris' inner rotation curve at z = for 
comparison. The star formation threshold has a significant effect on the mass distribution: a more prominent stellar bulge forms at early 
times in ErisLT and is responsible for the peaked rotation curve. 



found in lower-resolution simulations of Milky Way-sized 
galaxies that impose no restrictions on merger h istory: 
e.g., only one of th e eight galaxies simulated by |Scan- 
napieco et al. (2010| has a photometric D/T as large as 
U.t)8 (and six have D/T < 0.5), and only one ou t of the 
six ga laxies above A/vir = 10^^ Mq simulated by Brooks 



et al 
(TO 



2011) has a disk-to-total ratio comparable to Eris' 
which is offset, however, from the stellar mass- 



halo mass relation). 

3.3. Stellar Content 

Eris' total mass in baryons is Mi, = 9.5 x lO^'^ M©, 
corresponding to a mass fraction /f, = 0.12 that is 30% 
lower than the universal value (for the adopted cos- 
mology) of 0.175. Stars (and their remnants) comprise 
41% of all baryons within Rvir- the total stellar mass, 
M* = 3.9 X 10^° M0, is comparable to the valu e esti- 
mated for th e Milky Way, 4.9 - 5.5 x 10^° Mq, by |Flynn 
eraL]( |2006| . 

i'o make a bias-free comparison with the stellar mass- 
halo mass rel ation derived from the abundance matching 



technique by Behroozi et al. (2010) we adopt the follow- 
ing procedure! We ht the SUISy u, g, r, i, z broadband 
colors from the mock Sunrise im ages with the flexibl e 



stellar population synthesis code of Conroy et al. (2009): 
the fit assumes a Kroupa IMF and provides a ptioto- 
metric stellar mass estimate of A^* = 3.2 x 10^° M0 
(Conroy, private communication), 18% lower than the 
value directly measured in the simulation. The photo- 



consistently against the Behroozi et al 



age stellar mass-halo relation (which uses a U 



2010 


1 aver- 


a C 


habrier 



2003' IMF), free of IMF systematics, after offsetting aU 



j Behroozi et al.| (2010) stellar masses by 0.06 dex (to 
correct from Cnabrier to Kroupa IMF). The compar- 
ison, depicted in the right panel of Figure |4j demon- 
strates that Eris' implied "baryon conversion efficiency" , 
r/ = (A^*/Mvir) X (r^M/f^h) = 23%, is in excellent agree- 
ment with that predicted by the abundance matching 
technique. This contrasts with the recent analysis of 
many hydrodyna mic simulations of galaxy formation by 



Guo et al. (2010), who show that the great majority of 
them lock too many baryons into stars to be viable mod- 
els for the bulk of the observed galaxy population. Note 
that the intrinsic scatter in the stellar mass at a given 
halo mass is estimated to b e 0.17 dex, independent of 
halo mass ( [Yang et al.|[2009| . 

With a circular velocity at the radius, i?80 = 6.8 kpc, 
containing 80% of the i-band flux of Vgo = 210 kms~^, 
our galaxy lies close to the TuUy-Fisher relation of the 



Pizagno e_t_al. (20071 galaxy sa mple (see the le ft panel 
ot i-'ig. As discussed in |Pizagno et al.| (|2007), 

the TullyTisher uses Vgo as the primary velocity mea- 
sure rather than V2.2, the circular velocity at 2.2 disk 
scale lengths, since the former is less sensitive to the 
degeneracies of bulge-disk decomposition. The ratio 
^2.2/^200=214 kms~Vl29 kms-i = 1.66 in Eris, where 
V200 is the circular velocity at the radius enclosing a mean 
overdensity of 200pcrit (^200 = 177 kpc), is equal to the 
value s uggested by the dyn amical model for the Milky 



Way of Klypin et al. ( 2002 ) . It is also consistent with 



the recent measurements ot the v irial mass o f the Milky 
Way byiSmith et al.i 



plying 1^2.2/1^2 



200 



h ( |2007|) and |Xue et al.| ( [2008) , im 
1.48™ and ^2.2/^200 = l-67t°ii. 
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Fig. 2.- — Lef t panel: The optical/UV stellar properties of Eris at 2 = 0. The images, created with the radiative transfer code Sunrise 
) Jonsson|2006^ , show an i, V, and FUV stellar composite of the simulated galaxy seen face-on and edge-on. A Kroupa IMF was assumed. 
. iigtit panel: Projected face-on and edge-on surface density maps of Eris's neutral gas at 2 = 0. The color bar shows the neutral gas 
fraction. 



respectively]^ Note that while there is an unsettled dis- 
agreement among estimates ol the Galaxy's virial masf[^ 
Eris is likely 20% less massive than the Milky Way anc 
therefore the agreement between simulated and observed 
kinematic data should be confirmed in future simula- 
tions of slightly more massive halos and different accre- 
tion histories. Like the Milky Way, however, Eris is off- 
set relative to determinations using vari ous dark halo 
mass t racers for late- type disk galaxies by 'Button et al. 
(2010), who predict for typical dark matter halos with 
1^2.2 = 220 kms^i the ratio 1/^2.2 /V'2( 



l.llt°-J§ (2a) 



3.4. Gas Content 



* The V2 2/V200 ratios from | Smith et all |l2007| and |Xue et"aL] 
were computed by |Dutton et aI.] ( |^UlOp Irorn these data sets 
after converting different virial mass dehriitioiis and for an assumed 
Milky Way's V2.2 = 220 kms"!. 

^ Recent estimates of the virial mass of the Milky Way range 
from Mvir = l-0~']o'2 ^ lO'^^ Mq from blue horizontal branch kine- 
~^ ,-1-0.7 



matics i Xue et al.|2008 l to 1.2tol x 10^^ f^om studies based 



on the prop erties of the Large and Small Magellanic Clouds ( |Busha| 
|et al.|2011[ |. 



Eris' H I mass is A/hi — 1-9 x 10^ M©, comparable 
to the H I mass estimate d for the Milky Way disk by 



Nakanishi & Sofue ( 2003), but small e r than the value of 
^ 5 X 10^ Mq given by fWolfire et al] ( |2003[ ). The H i-to- 
stellar mass ratio, 1.9 x 10^ Mo/3.9 x 10^'^ Mq = 0.049, 
is equal to the median value observed in the GASS survey 
( Catinella et al. |2010 1 for galaxies of comparable stellar 
mass. Eris' H i disk extends out to about 15 kpc (6 stel- 
lar disk scale lengt hs), similar to the size of the H i disk 
of the Milky Way ( [Nakanishi fc Sofue||2003[ ). Clustered 
SN explosions create a large number or holes in the face- 
on H I distribution of Eris (Fig. [2]) due to bubbles of hot 
gas expanding perpendicular to the disk. These holes are 
mostly located within the bright optical disk and prefer- 
entially in regions of high star formation: they are kpc 
in size, as observed, e.g. , in the nearby low-in clination 
spiral galaxy NGC 6946 ( |Boomsma et al.||2008[ ). 

About 6.7 X 10^ Mq are found in a cold phase below 
T = 3 X 10'* K. This is comparable to the total mass 
of the atomic an d warm ionized medium inferred for the 
Milky Way (e.g. Ferriere 2001). The gas mass that is 
hot (T > 3 X 10^ K) and thus potentially X-ray luminous 
is Mx = 3.6 X 10*" Mq, 63% of the total gas content. 
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Fig. 3. — The ID j-band radial surface brightness profile of Eris 
at 2 = 0. This is well fitted by a Sersic bulge with index rts = 1.4, 
an exponential disk with scale length ij^ = 2.5 kpc, and a bulge-to- 
disk ratio B/ D = 0.35. The dust reddened, face-on 2D light distri- 
butio n created by Sunrise was analyzed with Galfit l |Peng et al. 
2002| follow ing a procedure similar to that detailed m | Wemzirl 
et al.| | |2009 jl . The "downbending" in the brightness exponential 
pronie at about 5 disk scale length and the surface brightness where 
the break occurs, 2 3.5 i-mag arcs ec^'^. arc characteristic of late- 
type spiral galaxies jPohlen &: Trujillo|2006[ l. 



For comparison, 12% of the gas is in the cold phase and 
25% is in a warm phase with 3 x 10"^ K < T < 3 x 10^ 
K. The fractions of cold, warm, and hot gas within its 
inner 20 kpc are 83.5%, 1.5%, and 15%, respectively. 
Hot gas within 20 kpc contains significant amount of an- 
gular momentum and is co-rotating with the cold disk. 
The hot gas baryon fraction, fx = MxlM^-„ = 0.046, is 
3.8 times smaller than the cosmological baryon fraction. 
This implies an average density for the hot gas that is 
3.8 times smaller than assumed in the standard "cool- 



ing flow" halo model a la [White fc Frenk (1991), and 
yields a factor of 14.5 smaller X-ray emission measure. 
Contrary to the standard assumption that hot gas fol- 
lows the radial distribution of the dark matter, the den- 
sity distribution of hot gas in Eris follows a "flattened" 



Px{r) 



-1.13 



power-law profile out to 100 kpc (see 



Fig. K^. This gives origin to an X-ray surface brightness 
profile that is not as sharply peaked as expected for hot 
halos with NFW profi les and tha t satisfi es the ob serva- 
tional constraints (e.g. Anderson & Brcgman p010| . The 
fiattened density profile we find is consistent with the 
resu lts of much lower-resolution simulations by |Crain et| 
al. (20101, who identified the reason for the more ex- 
tended gas distribution and weaker X-ray coronae in the 
entropy injection by SNe at z ~ 1 — 3. Only 10% of 
the gas in Eris is in a very hot phase above T = 10^ K: 
the mean density of million degree gas at i? > 70 kpc 
is n < 6 X 10~^ atoms cm"-^, which is well within the 
observa tional constraints from O vi absorption measure- 



ments (Sembach et al. 2003) in the halo of the Milky 
Way, and high enough to produce significa nt ram pres- 
sure stripping of dwarf spheroidal satellites ( Mayer et al. 



The observed dispersion measure (DM) of pulsars in 
the Large Magellanic Cloud (LMC) provi des another 
constraint to th e hot halo of the Milky Way ( Anderson fc 
Bregman 2010). Of the 11 pulsars discovered by 'Manch 
ester et al. (2006] ) in the direction of the LMC, 3 have di& 
persion measures below 45 cm~'^ pc and may be located 
within the Galaxy at a random position along the line 
of sight to the Cloud. The other 8 have dispersion mea- 
sures in the range between 65 and 130 cm^^ pc and are 
thought to be associated with the LMC, with the lowest 
dispersion values belonging to pulsars located on the near 
side of the Cloud, some 50 kpc away. This leads to an 
estimate for the dispersion measure introduced by Galac- 
tic free electrons towards the LMC of DM w 70 cm~'^ pc 



( Anderson fc Bregman][2bl0 1. To compare these values 
with prediction from the Eris simulation, we have cal- 
culated the mean integrated column density of free elec- 
trons between randomly positioned "observers" on a cir- 
cle in the stellar disk at a galactocentric distance of 8 
kpc, and a "pulsar" 50 kpc away. The mock pulsar was 
located at galactic coordinates I = 280° and h = —30° 
like the LMC. The predicted dispersion measure. 



DM = 



50 kpc 



,(/)d/ = (62±3) cm"3 pc (2) 



where is the electron density along the line of sight, 
is perfectly consistent with the data. Collectively, the 
above arguments indicate that the reservoir of hot gas 
around the Eris simulated galaxy appears to satisfy pul- 
sar dispersion-measure as well as X-ray surface bright- 
ness and O VI absorption constraints. 

3.5. Star Formation and Kinematic Decomposition 

The top left panel of Figure [6] shows the star formation 
history of all star particles identified within Eris' virial 
radius at z = (regardless of whether they formed within 
the main host or in satellites), and of its kinematically- 
decomposed present-day disk and spheroid. The de com- 
position technique follows Scannapieco et al. ( 2009 1 , and 
is based on the distribution ot orbital circularity pa- 



rameters. 



2z 
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of the simulated stars introduced by 
Here, jz is the angular momen- 



_ ( |2003[ ). 

tum of each star m the z-direction (i.e. the direction 
defined by the total angular momentum of all gas parti- 
cles within 5 kpc from the host center) and jc is the an- 
gular momentum of a circular orbit at the same radius. 
Spheroidal stars are defined as those that are not part 
of the rotationally-supported disk and therefore typically 
include bulge and stellar halo stars, as well as stellar bars 
if they are present. The distribution of circularity param- 
eters, shown in the top right panel of Figure |6j is charac- 
terized by two peaks: one at jz/jc — 1 that is indicative 
of the presence of a dominant cold disk in rotational sup- 
port, and a second one at jz/jc — corresponding to a 
modest hot spheroidal component dominated by velocity 
dispersion. 

The vast majority of disk stars in the Milky Way reside 
in a thin disk comp onent with expon ential scale height 



hz = 300 ± 60 pc ( Juric et al.||2008[ ). A study of the 
vertical structure or edge-on spiral galaxies finds that 
the scale height of their thin disk increases systematically 
with circular velocity as zq = 610 (14/100 kms~^)°-^ pc 
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log [km s-'] log M,., [M,^] 

Fig. 4. — Left panel: The i-band TuUy-Fishcr relation for the |Pizagno et ar] | |2007| l galaxy sample (empty squares with error bar s). Filled 
circle : The Eris simulation. Here Vgo denotes the circular velocity at the radius containing 80% of the j-band flux, as defined by lP izagno| 
|et al.| | [2007 i). Right panel: The stellar mass - halo mass relation at 2 = 0.1 from Bcliroozi et al. (2010), modified for a Kroupa IMF [empty 
squares with error bars). Errors bars include only statistical uncertainties. Filled circle: Ihe Ens simulation with a photometric stellar 
mass of M, = 3.2 X 10^^° Mq and a virial mass of Mvir = 7.9 X 10^^ Mq (see text for details). 
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Fig. 5. — The average dark matter {blue empty dots) and hot 
(T > 3 X 10^ K) gas {red empty dots) density profiles of Eris at 
z = 0. The solid lines show the best-fit NEW profile for the 
dark matter {upper curve) and the best-fit power-law profile (with 
slope —1.13) for the hot gas {lower curve). The best-fit NEW 
profile is characterized by a large halo concentration parameter 
c = iJvir/Rs = 22 as the dark matter halo contracts in response to 
the condensation of baryons in its center. 



( Yoachim & Dalcanton 2006 1, where zq is the scale height 
of a sech^ profile, zq ~ 2/1^ at large heights above the disk 
plane. Eris' kinematically-decomposed stellar disk agrees 
well with the above scaling: by fitting an exponential 
(sech^) profile to the simulation data, we derive a scale 



height of hz = 490 pc {zq — 860 pc) at a galactocentric 
distance of 8 kpc. 

Today, Eris is forming stars at a rate of SFR — 
1.1 MQyr~^, comparable to the value of SFR = 0.68 — 
1.45 Mf7^ yr"-'^, recently inferr ed for the Milky Way by 



Robitaille fc Whitney ( |2010 ) using Spitzer data. The 
star formation rate declines rapidly with redshift from 
a plateau value of ~ 10 MQyr"^ maintained between 
z = 2 and z = 5. SN feedback and photoheating by the 
ultraviolet radiation background efficiently quench star 
formation at z > 5. The rate of formation of spheroidal 
stars fades rapidly after redshift 3, while disk stars nearly 
triple their mass from z = 2 to the present. 

The star formation rate surface densities Ssfr and 
H I gas surface densities Shi (Kennicutt-Schmidt law) of 
Eris' disk are shown in the bottom-left panel of Figure [6j 
The simulation data were averaged over square patches 
750 pc on the side for comparison with the pixel-to- 
pixel observations at 750 p c resolution of 7 siJ iral galaxies 
from the THINGS survey fBigie l et al.][2008| . The same 
THINGS data are also plotted against the total gas sur- 
face density Sgas = Shi + Sh2 • The observed relation- 
ship between Sgas and Ssfr varies dramatically among 
and within spiral galaxies, and most galaxies show lit- 
tle or no correlation between Shi and Ssfr (Bigiel et 



al 



)8|. 



20081). Rather, there is a clear correlation between 
Ssfr (molecular Schmidt law), and gas at den- 
10 MfT, pc~^ is observed to be fully 



ani 

sities in excess of ~ iU IVI© pc 
molecular. At solar metallicities, this corresponds to the 
column of atomic hydrogen required to shield a molecu- 
lar region against p hotodissociation ( Krumholz, McKee 
I& Tumlinsonl'2009V While the total mass of cold gas in 
Eris is comparable to the sum of the atomic and molec- 
ular gas mass of Sb galaxies such as the Milky Way, we 
do not model directly the formation of H2 molecules. In- 
deed, as clearly shown in Figure |6j even at Eris's high 
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Fig. 6. — Top left: Star formation history of all star particles identified within Eris's virial radius today. Black filled dots: total 
star formation rate {top panel) and stellar mass {bottom panel) as a function of redshift. Blue filled dots: same for disk star particles 
identified at ^ = 0. Red filled dots: same for spheroid star particles identified at 2 = 0. See the text for details of the disk-spheroid 
kinematic decomposition. Top right: Stellar mass fraction as a function of the "orbital circularity parameter" jz/jc, describing the degree 
of rotational support of a given stellar particle, for Eris at z = {solid line), Eris at z = I {short-dashed blue line), and ErisLT at z = 1 
{long-dashed red line). The prevalence of stars in a centrifugally-supported thin disk manifests itself in a sharply peaked distribution about 
unity. Bottom left: i^sFR versus Ejji for Eris' disk at z = {black filled dots). The simulation data were averaged over square patches 750 
pc on the side: some discreteness effects associated with the limited resolution of the star formation timescale can be seen at low values of 
^SFR- Every dot represents one sampling point. Note that our simulations do not model the formation of molecular hydrogen. The blue 
emvtv d ots show the pixel-by-pixel Sgpji data as a function of Shi (a-t 750 pc resolution) for 7 spiral galaxies from the THINGS survey 
( [Bigiel e t al. 2008). The same THINGS data are plotted against the total gas surface density Sgas = Shi-I-Sh2 {gray empty dots). Bottom 
right: Evolution of the baryon fraction within the virial radius for Eris {blue filled dots) and ErisLT {red filled sguares). In the adopted 
cosmology, the cosmic baryon fraction is = 0.175. 
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resolution we do not capture kpc-sized regions with gas 
surface densities in excess of the characteristic shielding 
column. Furthermore, our gravitational softening is still 
large compared to the size of typical molecular clouds, 
which are a few tens of pc in size, and our gas density 
threshold for star formation is still somewhat below the 
density of real molecular clouds. Both effects contribute 
to further smooth out the high density tail of the total 
gas distribution. Yet the figure shows that our simulated 
disk is forming stars in the same range of Ssfr " Shi 
values observed in spiral galaxies. 

The global atomic gas depletion timescale, i.e. the time 
needed for the present rate of star formation to consume 
the existing atomic gas reservoir, is tni = Mhi/SFR = 
1.9 X 10^ Mq/I.I M0yr-i = 1.7 Gyr in Eris today. As 
this is significantly smaller than the Hubble time, star 
formation rates and gas fractions must be set by the bal- 
ance between gas accretion from the halo and stellar feed- 
back. The depletion time for atomic gas ob served in the 
COLD GASS survey ( |Saintonge et al.|2011[ ) IS on average 
3 Gyr, with a large scatter from one galaxy to another. 

4. DISCUSSION 

Approximately 70% of bright spirals have B/T < 0.2 
( Weinzirl et al.|2009 ). Recent attempts to generate such 
disk-dommatecf galaxies in cosmological simulations have 
failed to reproduce simultaneously their observe d mor- 
phologies as well as their baryonic /stell ar content ( Scan- 
napieco et al. 2010 Agertz et al.||201l" |. Our very high- 
resolution Eris simulation appears to form a close ana- 
log of a Milky Way disk galaxy by capturing a realistic 
inhomogeneous ISM in which star formation occurs in 
high density regions of mass comparable to that of giant 
cloud complexes. Contrary to the "inefficient star forma- 
tion" prescription adopted by Agertz et al. (2011 ), this is 
achieved with a strong localized SIN feedback and a high 
(10%) Schmidt-law eflSciency. We stress that our effi- 
ciency parameter is simply phenomenological, and does 
not have a direct relation to the true star formation ef- 
ficiency within giant molecular clouds, which is the end 
result of all the processes regulating star formation in- 
cluding feedback, and concerns scales (tens of parsecs) 
that are not yet resolved in cosmological simulations. 
While the chosen star formation and feedback prescrip- 
tions combine to expel a significant amount of baryons 
from the system, they preferentially remove low angular 
momentum material and leave plenty of cold gas avail - 
able for disk star formation (cf . [Scannapieco et al.|2009[ ), 
thus allowing a good fit to the I'uiiy-l'isher relation (cf. 
Piontek fc Steinmetz|[2011[ ). 

It is interesting to compare the properties of Eris and 
ErisLT at redshift 1 (see Table 1): 

• The ErisLT control run produces a galaxy resem- 
bling an early- type Sa spiral with B/D = 0.42 
(cf. Eris' B/D = 0.31), closer to (but stih on 
the low side of) the typical outcome of previously 
published cosmological simulations of disk forma- 
tion ( e.g. |Governato et al.^2010, Scannapieco et al. 
2010). Its rotation curve peaks at 308 kms~^ (cf. 



Eris and close to the universal value. The differ- 
ence in the baryon content of the two simulations 
is established at very early times (see the bottom 
right panel of Fig. [6]). The baryon content of 
ErisLT increases from z — 10 to z = 5, as its 
dark matter halo grows from Afvir = 3 x 10^ Mq to 
Mvir = 6 X W^Mq. At higher redshifts (smaller 
progenitor mass), the collapse of baryons is heav- 
ily suppressed by the ultraviolet radiation back- 
ground. Eris, by contrast, maintains a relatively 
low baryon fraction, between 64% and 73% of the 
cosmic value at all redshift z < 10. Note that these 
values are higher than the value measured at the 
present epoch in the bulgeless d warf galaxy simula- 
tion of Governato et al. (2010), which is only 30% 
of the cosmic traction. 

• While ErisLT was run with a star formation ef- 
ficiency that was half of that adopted for Eris, 
its stellar content at z = 1 is 20% higher, show- 
ing that the most important parameter in the star 
formation recipe is not the efficiency but rather 
the star formation density threshold. Indeed, by 
allowing the gas to reach higher densities before 
turning into stars, the ISM develops a more inho- 
mogeneous structure, with important consequences 
on the large-scale pattern of star formation in the 
galaxy, and, as a by product, on the effect o f super- 
novae feedback (e.g. Governato et al.||2010 |. 



Eris' 237 kms ^) and declines steeply within few 
kpc from the center (Fig. [T]). 

The baryon fraction in ErisLT is higher than in 



At z = 1, about 14% of the total gas content in 
ErisLT is in a cold phase, 19% is warm, and 67% is 
in a hot component. The same fractions are 30%, 
18%, and 52% in Eris, i.e. there is 2 times more 
cold gas in Eris than in ErisLT. The difference in 
the density distribution and thermodynamic state 
of the ISM in the two simulations is clearly seen 
in Figure [7) which shows the local gas density ver- 
sus radius of each gas particle at different epochs. 
The horizontal green line in each panel marks the 
adopted star formation threshold, and particles are 
color coded according to their temperature. Stars 
are born only in the blue gas particles above the 
green lines. Because of the increased threshold, 
fewer regions are forming stars at any given time 
in Eris. According to equation ([T]), however, within 
these cold dense clumps the number of massive 
stars per unit gas mass born at threshold scales as 
Y^nsF, i.e. it is a factor (5/0. 1)^^^ = 7 times larger 
in Eris versus ErisLT. As these stars explode as 
SNe, the energy injected at high redshift in Eris's 
ISM is enough to disrupt and unbind the star form- 
ing regions, generating strong outflows that leave 
the galaxy and reduce its baryonic content. 

Figure [7] also shows that, at z = 5, star formation 
can only occur in Eris' inner densest regions within 
10 kpc from the center. By contrast, in ErisLT, 
star formation is spread out over a significant frac- 
tion of the host galaxy, extending well beyond 20 
kpc. At z = 3, when outflows are stronger due 
the high star formation rates triggered by merg- 
ers, it is the cold, low angular momentum gas at 
the center that is preferentially removed in Eris, as 



10 



The Eris Simulation 




5 10 15 20 5 10 15 20 

t [kpc] t [kpc] 





5 10 15 20 5 10 15 20 

t [kpc] t [kpc] 



Fig. 7. — Properties of the gas distribution for different star formation thresholds and different redshifts. The local gas density measured 
around each SPH particle is plotted as a function of its radial distance from the galaxy center for Eris {lejt panels) and ErisLT (right panels) 
at 2; = 1, 3, and 5 (from bottom to top). Horizontal green lines mark the minimum gas density for star formation in each run. The color 
coding highlights cold (T < 3 X lO" K, blue), warm (3 X 10* K < T < 3 X 10^ K, yellow), and hot (T > 3 X 10^ K, red)) gas particles. Star 
formation only occurs in blue gas particles above the horizontal green lines. 
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previously fo und for dwarf galaxies by |Governato| 



et al.| (2010). Such outflows are weaker in the low 
threshold simulation, where star formation is more 
diffuse and the energy injection from SNe is more 
evenly deposited in the ISM: cold gas accumulates 
in the inner regions already at very high redshifts 
and continues to turn into stars unabated by feed- 



back (see als o Ceyerino fc Klypin| 2009 Robertson 
&: Kravtsov||2008|). At lower redshitt ErisL'i' has 



now consumed its low-angular momentum cold gas 
in star formation, while Eris has preserved a higher 
cold gas fraction due to the more effective regula- 
tion of star formation via feedback. By z = 1 both 
galaxies have formed roughly the same amount of 
stars, but 90% of the gas within 20 kpc from the 
center is still cold in Eris, compared to only 57% 
in ErisLT. The distributions of stellar angular mo- 
menta are also different in the two cases (see Fig. 
[6]). Gas being blown out at high redshift has a 
systematically lower angular moment um than the 
gas t hat gets accreted at later times (Brook et al. 



2010), and this explains the flatter rotation curve 
ot Eris relative to ErisLT. 



It is fair to point out that, while the success of our 
Eris simulation in matching the observations appears to 
be linked to the ability of correctly following star forma- 
tion in an inhomogenous ISM and regulating it with SN 
feedback, many avenues remain unexplored and require 
further investigation. Simulations of even higher resolu- 
tion, approaching the true gas densities reached in star 
forming giant molecular cloud complexes (about a fac- 
tor of 10 higher than adopted here), are needed to test 
the convergence and robustness of our results, and are 
in the making. The density threshol d should depend on 
metallicity and therefore on redshift ( Gnedin et al. 2009 ) , 
which may affect the structure of the ISM and SJN teed- 
back differently in progenitors at different epochs. The 
feedback model that we use is still phenomenological, 
and the actual mechanism of outflow generation may re- 
quire the combination of mor e than one effect to suppor t 
a large-scale blastwave (e.g. Ceverino & Klypin 2009). 
Lastly, the simulations presented in this paper neglect 
cooling by metal lines at temperatures above lO'' K. At 
the mass scale considered here, cold flows are mostly re- 
sponsible for the assembly of the star forming disk even 
at low redshift, as opposed to the cooling flow of the 
hot halo mode (Brooks et al. 20091. This suggests that 
the details of the cooling function for gas above 10^ K, 



namely in the hot mode, are not important fo r the as- 
sembly of the disk. Piontek & Steinmetz (20111) find that 
metal cooling gives origin to a stronger burst of star for- 
mation at high redshift in the progenitors of Milky- Way 
sized galaxies and thus to bigger bulges. This result, 
however, was obtained with the canonical low star for- 
mation density threshold (nsF = 0.1 atoms cm~^), and it 
is unclear in which temperature range the effect of metal 
cooling is most crucial. These authors also show that 
the augmented bulge can be suppressed by boosting the 
effect of (thermal and kinetic) SN feedback. It is con- 
ceivable that, by further increasing the density threshold 
for star formation toward actual giant molecular cloud 
densities, along with modeling the formation of molecu- 
lar hydrogen, the ISM may become increasingly clumpy 



and dense, populating the high surface density tail of 
the Kennicutt-Schmidt relation that remains currently 
unresolved in Eris. As a result, one may expect that 
heating and outflows from even more localized SN feed- 
back may become stronger in such high density regions, 
and eventuall y offset the impact of i ncreas ed cooling, as 
suggested by Piontek & Steinmetz (2011). We plan to 
explore these issues in the next generation of simulations 
with increased resolution, a higher star formation density 
threshold, and the inclusion of molecular phase physics. 
In addition to the effect on disk and bulge, metal cooling 
at r > lO"* K may also have an impact o n the den- 
sity profile and t emperature of the hot halo. [Piontek 
Steinmetz (2011) find a significant reduction of the mass 
ot the hot halo when metal cooling is included, under 
the assumption of a gas metallicity of Z = 0.5Zq. This 
is higher than the metallicity measured in spiral galax- 
ies with extended X-ray and Ha emission a round their 
disk, which is in the range Z = 0.01 — O.IZq (Rasmussen 
et al. 2009), suggesting that the effect may have been 



overestimated. The mean metallicity of hot gas in Eris 
closely matches the latter observati ons, b eing o n aver- 
age Z = O.OSZq (with Zq = 0.0194, |A^rs fc Grevesse 
1989|. 



The last major merger in our simulations occurs at 
2 ~ 3, and therefore Eris is expected to show some offset 
from the observed structural parameters of the average 
spiral galaxy. The same likely applies to the Milky Way 
itself, which indeed closely resembles Eris. Whether or 
not the good match with the properties of typical spiral 
galaxies is related to the fact that we have selected a 
particularly quiet merging history, in which outflows shut 
off early as the rate of star formation drops after the last 
major merger, will have to be investigated. 

At z < 3, the evolution of the B/D ratio is non- 
monotonic and it is seen to decrease following a major 
merger at z = 3 and a minor merger at z = 1, and to 
grow secularly at lower redshifts. The details of this evo- 
lution will be the subject of a forthcoming paper (Guedes 
et al. 2011c, in preparation). More typical galaxies un- 
dergoing significant mergers at later times may develop 
smaller bulges due to the more prolonged effect of su- 
pernovae outflows, which could explain why 11 out of 19 
nearb y massive spirals show no evidence for a classical 
bulge ( Kormendy et al.||2010 l). If true, the trend would 
be at odds with the standard picture in which mergers 
lead to earlier type objects. This would be an important 
new ingredient, perhaps complementary to the finding 
that gas-rich mergers assist the growth of larger disks 
( [Hopkins et al.||2009[ [Governato et al]|2009[ ). 
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